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ABSTRACT 


Large orbiting space structures are expected to experience mechanical vibrations arising from several 
disturbing forces such as those induced by shuttle takeoff or docking and crew movements. In this paper, we 
consider the problem of modeling and control of large space structures subject to these and other disturbing 
forces. The system consists of a (rigid) massive body, which may play the role of experimental modules 
located at the center of the space station and flexible configurations, consisting of several beams, forming the 
space structure. A complete dynamic model of the system has been developed using Hamilton s principle. 
This model consists of radial equations describing the translational motion of the central body, rotational 
equations describing the attitude motions of the body and several beam equations governing the vibration of 
the flexible members (platform) including appropriate boundary conditions. In summary the dynamics of the 
space structure is governed by a complex system of interconnected partial and ordinary differential equations. 

Using Lyapunov’s approach the asymptotic stability of the space structure is investigated. For asymptotic 
stability of the rest state (nominal trajectory) we have suggested feedback controls. In our investigation 
stability of the slewing maneuvers is also considered. 

Several numerical results are presented for illustration of the impact of coupling and the effectiveness 
of the stabilizing controls. This study is expected to provide some insight into the complexity of modeling, 
analysis and stabilization of actual space structures. 


1 INTRODUCTION 

Structural vibrations in future large space systems(Fig.l) such as space antennas, space platforms, space 
stations, or deployed flexible payloads attached to the space shuttle orbiter, and their interaction with the 
other members of the system have become a major concern in design of reliable systems satisfying stringent 
stability requirements. During the past few years, considerable attention have been focused on the development 
of mathematical model and stabilizing controls for such systems[l-13]. 

The most natural model for flexible space structures is given by a hybrid system of equations which is a 
combination of ordinary differential equations for rigid parts and partial differential differential equations for 
flexible members. Hybrid models for some simple structures have been considered m [5-11J. Recently Limll 2j 
developed a complete dynamic model, which includes the orbital dynamics, attitude dynamics and equations 
for vibrations of flexible members and all the relevant boundary conditions, for flexible space structures. Here, 
in this paper, based on the dynamic model we develop a control scheme to suppress the vibration induced by 
slew maneuvers in space stations. 

The paper is organized as follows. In the next section we present the equations of motion for large 
spacecraft derived in [12]. In section 3, based on Lyapunov’s approach, asymptotic stability of slew maneuvers 
for the system using feedback control is considered. For illustration, in section 4, we present some numerical 
results demonstrating the effectiveness of the stabilizing controls for vibration suppression associated with slew 
maneuvers. Finally, the concluding remarks are presented. 
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Fig.l The U.S. Space Station. 
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2 DYNAMIC MODEL OF FLEXIBLE SPA CE STATIONS 


2.1 Introduction 

We describe the dynamics of the space station by 
Three Types of Motio n: 

- Rigid body translation perturbing the orbit. 

- Rigid body rotation perturbing the orientation of the system. 

- Vibration of the elastic members causing elastic deformation of the system. 

To derive the dynamics we use the following coordinate systems: 

2.2 Reference Coordinate Systems 

- Body Coordinate System Sb ' (h,jb,kb) 

- Orbital Reference Coordinate System 5 r : (*rj jn*r) 

- Inertial Coordinate System 5/: 


Angular Velocities 

Wb = (t v Xj w y ,w z y: angular velocity of body frame w.r.t. the S r . 
w r — {wxi wy , tn^)':angular velocity of S r w.r.t. 5/. 
w = W b + W r = (w\,W 2 , u? 3 ) / :angular velocity of 5 b w.r.t. 5/. 



n = R + fi , with r,- = ( Xi,yi , z,*)' = D io + Di , 
where R = ( X,Y,Z D, 0 = (i <0 , y< 0> and Di = Q, ifc, *<)'. 

= (0,/,), i=l,2,...,6 for /.(the length of the beam i). 

D\ 0 — (si, + £i, «i v , 0/, Di, — («2* +^2, *2„>0) , 

Dz o = ( — ($3* + £3)1*3, >0/, = ( — ( s 4* + £4)1 S 4,1 0) , 

Z? 5 „ = (® 5 , i — s 5, +^5>0) , ) De, = (*6, *6, + £e>0) > 


( 1 ) 


( 2 ) 


(3) 

(4) 

(5) 


Jt 
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Fig.2 The space station and reference coordinate systems. 
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2.3 Inertia Tensor 

The beam inertia tensor is given by 


where 


*-xy 


l V* 


/ Ixy *xz \ 

1 ~Jxy lyy ~Iyz 1 
\ —Izz lyz Izz ) 

(6) 

E / 

(7) 

E / «(*? + *?)«%. 

(8) 

E/ 

(9) 

E / 

(10) 

E / PiZiZidti, 

TZ tJfii 

(11) 

E / PiViZ&i- 

(12) 


2,4 Dynamic Mod el of the Space StatiQn 

Let the Lagrangian of the whole system be denoted by L and external work by W . Then using the 
extended Hamilton’s principle we obtain the relation[12,13]: 

6 [\L+W)dt 

Jt i 

(6Di)oI a d$i+l 4 ^dt 


= 0, 


where 




d? R d f ( dR . , n \ j , Gm e m T R 

+ sgy„. \ it + - “ x r * + nr ™ 1 + ~w~ 

l! = T - !{<'• °“> + x U? x * + A ) x (1« 

6 ' 

■ s{K£*)4H«, - J-'- ■ w 6 


(13) 


(14) 


(15) 

(16) 


(17) 
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Let 


/ EI it 0 0 \ 

EIi=\ 0 EI iy 0 , i=l,2,.. .,6. 

\ 0 0 El it ) 

F = f (0 ,F t „F t 'Y, for *=1, 2,3, 4, 

‘•"l (F it ,0,Fi.y, for »=5,6. 

For boundary conditions we introduce the following assumptions: 

Assumptions for the Joints and Beams: 

(i) The structures are in xy plane in undisturbed state, 

(ii) The beams 1,2, 3, 4 are clamped at the joints J \ , J 2 

(iii) Displacements of the beams 1,2 ,. ..,6 are small, 

(iv) The beams 5,6 are rigidly jointed with the other four beams as shown in Fig. 2. 


(18) 

(19) 


For convenience of presentation we use the following notations: 

ati = <1 j x ri + x Di + u> x (u> x r,), i = 1,2, 6. 


Fi= [ ( Fi ,,/U'i 
4 ‘~l (Fi.'FJ, 

(° 1 

Vo 0 1 )' 


Qi = { 


for i=l, 2, 3,4, 
for i=5,6. 

for i=l,2,3,4, 


( 001)’ for i=5,6. 


mass of beams : 774 


-U 


dm,-, mass of the body = m 9 


mass density of the beam : pi = />*(&), m r = m s + 
6 




= 1 fi 

*=1 j 


x (c 4; x ri)drrii 


I T = I, + h f 1$ : inertia of the body h : inertia of the beams. 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


Since the variations 6R, 60 and 6D{ are arbitrary in (13), I; , : = 1,2,3 are all zero provided J 4 = 0 from 
the boundary conditions. Hence we obtain from this fact the following equations of motion[12,13]: 


Ti=0, that is, 


d 2 R ■sr~' f ^ ■ 1 ,. Gm e m T R 

mT ~dtF + x r *) +w x x r ‘) + +2w x D ijdt,i + — — = F, 


(26) 


J 2 = 0 , that is, 


I T o u + ui x (I T o u>) + i b o u> + J n Pi(r> x 

+ p '( f ' x Fi)d^ + YlJ a Pi x ( f< x A)}de, = T , 


( 27 ) 
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and for & € fii, £j € *=1 ,2,3,4, j— 5,6, 1z— 0, that is, 


a 2 

^0< 2 

(:) 

+ ( 

f d 2 R\ 

, a ’ + H 2 ) 

l+ l?( 


o > 

> 

\ — 

(“))= n " 

(28) 

a 2 

Pi a< 2 

(!) 

+ Pi Qj 

( cPR' 

P + dJ 

) + lj( 

( £ 5" 

0 

El iti 

) AC/ 


(29) 


Let Cg denote the coordinate transformation matrix from body to inertial frame and (a, , a ; , a^)' _ (Gg) 
£iR. Then the equations of motion given in (26-29) can be written in explicit form for computer simulation 
as follows. 

BODY DYNAMICS!" TRANSLATION) 


m. 


d 2 


' Y ' 






/M<)\ ( 

MO = 
WO/ \ 


F i' 

F 2 

f 3 


where 


dyi . 

■"•■ST* 


/i(o = ^2 [ pi \ ( 2,W2 - y,W3 ) + 2 ( W2 ir - 1 

i= l Jfl * *- 

d 2 Xi \ 

+ Wi w 2 y< + u>iu> 3 ii - (ti > 2 + «>!)*,• + -Q-j- 

® r ( . . dxi dzt , 

h{t) = Z2J a Pi\*i W 3 ~ * iWl + 2 ( W3 ~dt ~ Wl ~dt^ 

q 2 \ 

+ w 2 w 3 Zi + w 2 w x Xi - (tt >3 + u> 2 )yi + -Q-f- |d&, 

6 m S Q % ’ 

/s(0 = YL I £i ™ 2 + 2 ( u ' 1 ~it - W2 ~dt ^ 

d 2 Zi \ 

+ ufewiSj + tv 3 w 2 yi - (w? + w\)zi 4- 


l (x\ _ Gm e m T 

W) - 

» {i \ Gm t m r 

M0--| S j3“> r . 

. o\ _ Gm t m T 
A e(0- 1^3 


PODY DYNAMICSfROTATION) 


It | w’a 1 + 

1 W 3 , 


0 — w 3 

U?3 0 

— U>2 Wi 


W2 \ / Wl \ ( f 4 + fl\ { 

— Wl 1 It I W2 I + I /s + fa I = I 

0 / \ u 3 / V /6 + /9 / \ 


Tx 

T 2 

t 3 


where 


(30) 


(31) 


(32) 


(33) 

(34) 

(35) 

(36) 

(37) 


m = t j *{*f£ - * w + 2 »>*il + - J -»*7r - ^irK < 38 > 

i=l Jn < K 
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hit) 


/e(0 


-U .4 


_ d 2 Xi _ d 2 Zi _ dzi _ dxi _dt/i _ 

- W - Xi — + 2u, 2 z,— + 2 w 2Xi — - 2w 3 z i — - 2w lXi 


dt J^*’ 


. d 7 yi _ d 2 Xi _ dxi _ dyi _ dzi _ dzi . 

' W - yi ~w + 2w3Xi ~m + 2w3yi ir - 2wiXi ~m - 2w2Vi ~m >**• 


(39) 

(40) 


m = 

hit) = 
hit) = 



BEAM DYNAMICSfBEAM VIBRATION) 


„ 52 f* 

p, dt 2 v 

-2 Ul \d(yi 
:i)+ Pi \ 2 u>j 0 )dt\zi 

i\. d 2 ((EI iy 0 \ d 2 /y,\\ 

J EI iz ) dt 2 W J 


+ Pi ( 7 W ? " w2 - W1 + W2 " 3 Pi xi ( 0)3 + L W1W2 ) 

\Wl+W 2 W 3 -wf - U>$ J \Zi J \-U) 2 +W1U3 J 



for ti € * = 1.2, 3 , 4 , 

d 2 (x 

Pi dt 2 c 

■?)!(: 

■A + 0 L((e*j* 0 

j) dtj\\ 0 EI iM ) dt] \*i) ) 


( tv 2 -wl W 2 + W3W1 'N fxj \ f —W3 + W 2 Wi \ 

Pi \ — t b 2 + W3W1 -w% — w 2 J \Zj J Pi 2 \ thi + W 2 W3 ) 


+ P & C *W - {%]) ’ 

for tj € Hj , j = 5 , 6 . 


(41) 

(42) 

(43) 


(44) 


(45) 


BOUNDARY CONDITIONS FOR BEAM DYNAMICS 

Since the beams(l-4) are clamped to the central body and the end beams(5,6) are rigidly jointed to the 
first four beams as shown in Fig. 2, the following boundary conditions must hold. 

Clamped at £,* = 0, i= 1 ,2,3,4 for the joints J\, J 2 , J 3 , Ja- 


y,(0, <) = 0, *,(0,0 = 0, (46) 

|g-(°.0 = °. ||(<M) = 0, (47) 

Continuity of displacements at the joints 

*5(f5,0 = 0, * 1 (h,t) = zs(ls,t), at Ji, (48) 

Z5(0,t) = 0, z 2 (h,t) = *5(0, t), at J 6 , (49) 

*6(^6, 0 = 0 , Z 3 (l 3 ,t) = Z 6 (l6,t), at J 7 , (50) 

X 6 (0,0 = 0> z *(U,t) = *6(0,0. at J&y (51) 
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(52) 


Continuity of slopes at the joints 


Equilibrium equations: 


|g(M> = -fg(°, <),«<*. 

!^,o= g«..o 

% (,, ' <)= I; 10 ' 0, al Js ’ 


at Js, Eh z -7^(h,t) = 0 , 

a<J 6 , £/ 5 ^( 0 ,<) = 0 , EI 2 ,^(h,t) = 0 , 

a<J 7 , 2?J6 z^(/6,<) = 0, EI 3 ,j^{l 3 ,t) = 0, 

at J s, EI$ t 0£2 6 "(®’*) = EI^ t -^^{}A,t) = 0, 

Ehy^fVut) = 0 , EI 3y ^(h,t) = 0 , 

EI 2 y g^(l2,t) = 0, Elly ^3 (^4,<) = 0> 

EI u + EI *z^(ls,t) = 0 , 

£7 *W (,8,<) + Fl6i W (,6,t) = 0> 

£J 2 ^(/ 2 ,<) - = °. 

Eh,^(U,t) - Eh,^( 0 ,t) = 0 , 

- Ehr^rih,t) = 0 , 


Ehy-^rihJ) + Eh*-jg(k,t) = 0 , 
Ehy^(h,t) + Eh*^( 0 ,t) = 0 , 

(/4><) ” El6t ^[ ( - 0 ’ t) = °- 
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3. STABILIZATION OF SPACE STATIONS 


3.1 The Reference Trajectory 

Let Rc denote the reference trajectory(nominal trajectory) governed by the equation 


<PRo Gm^Ro 

df 2 |f?o| 3 

with the initial conditions: 

dRo ( _ dRo\ Gm e Ro 

T x r x T) = W 


(70) 


(71) 


Denoting the excursion of the radial vector R around the nominal trajectory by R _ R Ro — (^>^>^0 we 
obtain from eqs.(26) and (70) the perturbed radial dynamics: 


<PR 


dt 2 



Gm t m r R 


(72) 


3.2 The Rest State 

For stability of the system subject to external disturbances we consider the rest state: 

U>1 =0, tl> 2 = 0, u>3 = 0, 

X(t) = Y(t) = Z{t) = 0, Vy(t) = e 2 (<) = W3(<) = 0, 
and for fc = 5,6; j = 1, 2, 3,4; i = 1,2, ...,6, 

= *»(£«> 0 = 0 , ~ ~ = 


(73) 

(74) 


(75) 


where v = ^ = («i, v 2 ,tJ 3 )'. 


Theorem 1 (Distributed Control) 

Consider the perturbed system described by the equations (30-69) around the reference trajectory. Sup- 
pose that the controls applied to the system are given by the feedback law: 


T = (T r — k\ Wi, Ty — * ^2^2) Tf — ^ 3 ^ 3 ) » k$ > 0, 


Fb, —Qi ( 


4.. i, -d.,^V 


“ at 


dt' Ai ‘ dii dt)' 


4 1 .4„4,>0, * = 1,2, ...,6, 


F, = (-e^i, -e 2 t3 2 , - 63 ^ 3 )', 


ei,e 2) e 3 > 0, 


(76) 

(77) 

(78) 


where (A,,, A <J( , A*,)' = and (T X) T 9 ,T t )’ = -f^(E?=iJ n< *^.) x *o- Then the s y stem is 

asymptotically stable (in the sense of Lyapunov) with respect to the rest state (73-75). 

Proof see the refs.[12,13]. 


Remark 

Defining the state variable x appropriately and incorporating the boundary conditions in the differential 
operator one can rewrite the system equation as an ordinary differential equation in a Banach space (probably 
Orlicz space) as follows 

x = Az + F(x,x), (79) 

where A can be proved to be the infinitesimal generator of a C 0 - semigroup in the Banach space and F is a 
strongly nonlinear operator having polinomial growth. This is a descriptor system in an infinite dimensions 
space and very little is known about these systems concerning the existence and uniqueness of nonlinear 
semigroups. 
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4 NUMERICAL RESULTS 


For numerical simulation we assume that (i) the bus inertia tensor is diagonal, (ii) the flexible members 
consist of three beams (i.e., beams 1,2 and 5 in Fig. 4) and each beam is uniform. 

The following data and parameters were used. 

System Parameters : 

7* =diag( 7i, 7 2 , ^3) — diag(31750, 5000, 33450) slug ft 
the length of the beams: /1 =/ 2 =187.7 ft } 1 5 =66 

Flexural rigidity: EI$ X =EI$ Z =10 5 lb ft 2 , EIi V =EIi z =3565.5 lb ft 2 ) for i=l,2,5, 

Mass density: pi =8.25 x 1CT 2 slug/ ft, 

D lo =(3+ ft, 33, 0)', 77 2o =(3+ ft, -33,0)', D 5o =(190.7, -33+ ft, 0 )'ft } m s =300 slug. 


Initial Conditions : 

uq (0)=0.03, u> 2 (0) =0.02, u; 3 (0) = 0.01 rad/sec and for ft £[0,/*] and <=1,2, 5, ^(ft>0) =*^Ki)0) 
= ^(ft, 0) =0, x<(ft,0)= y,(ft,0) =*,(ft, 0) =4>o (ft), where <f> 0 satisfies the boundary conditions (equations 
3.97-3.120). 


Control Gains : 

=30000, k 2 =10000, k 3 =100 for t £[0,2.5] and k 3 =50000 for t €[2.5,10]; d h =d h =d< 3 =0.05 for i= 
1,2,5; ei=9000, e 2 =3000, e 3 =6000. F x (t) =F 2 (<) =^(0 =0. 


Slew maneuver using Bang -Pause -Bang control with the parameters: 
7^ = 10035, < #l =0.5sec, < #a =2.0sec, f #3 =2.5sec. 


slew control 


T.. 







i 

■>, t 

*3 

-r tl 








Fig.3 Bang- Pause-Bang Control. 

Detailed numerical results showing stabilization of the various state trajectories were obtained from the 
simulation and presented in Figs. 4-23. Figs .4 and 5 represent the slew angle and the angular rate of the 
spacecraft corresponding to the slew control (Fig.3). We note from the Fig. 4 that slew maneuver was suc- 
cessfully achieved by the slew control commands. As one can see from Eqs. (30-45), the body translation, 
body rotation and vibration of beams are strongly coupled. This implies that any perturbation in one part of 
the system will induce disturbances in the other members of the system. Hence the slew maneuver induced 
significant perturbations leading to beam vibrations and body rotation. This is clearly observed from the 
responses without controls. From Figs. 4-23 it is clear that without stabilizing controls beam vibrations, body 
oscillations and radial excursions persist or grow. However, with application of the proposed feedback controls, 
oscillations induced by the slew maneuver are effectively eliminated throughout the entire system. 
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Fig. 6 Bus angular velocity wj. 
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Fig. 7 Bus angular velocity wj. 
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TIME (SEC) 
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Fig. 10 Beam displacement jfa{Jj,0* 


Fig II Beam velocity 
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WITHOUT CONTROL 


WITH CONTROL 
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Fig. 10 Beam displacement 


TIME (SEC) 
Fig. 13 Beam velocity 


WITHOUT CONTROLS 


WITHOUT CONTROL- 


WITH CONTROL - 


WITH CONTROL 


TIME (SEC) 

Fig. )4 Beam displacement Xi(l^,t). 


TIME (SEC) 
Fig. 1 5 Beam velocity 
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Fig. 20 Radial perturbation k. 
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Fig /2 IUdi»l p«turb»tion Z. Fig 22 ReUti.e bemm energy. 
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■V CONCLUSIONS 


In this paper, suppression of vibration induced by slew maneuver in flexible space stations has been 
considered. Based on the dynamic model developed by the authors[13], asymptotic stability of the sys em 
subject to perturbation is investigated. The rotational perturbing forces were applied to the system and their 
corresponding stabilizations have been demonstrated. From the numerical results it is clearly observed that 
fi)if the disturbances following external perturbing forces persist, then in the absence of proper controls, these 
small motions may build up leading to instability of the entire system, (ii)during the slew maneuver, vibration 
in the beams and oscillations in the angular velocities of the body are induced and it has been shown that the 
stabilizing control can effectively eliminate the oscillations and stabilize the entire system. 
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